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Abstract 



Using Dirac's approach to constrained dynamics, the Hamiltonian formu- 
lation of regular higher order Lagrangians is developed. The conventional 
description of such systems due to Ostrogradsky is recovered. However, un- 
like the latter, the present analysis yields in a transparent manner the local 
structure of the associated phase space and its local sympletic geometry, and 
is of direct application to constrained higher order Lagrangian systems which 
are beyond the scope of Ostrogradsky's approach. 



1 Introduction 



Canonical and path integral quantisations of systems whose dynamics is de- 
scribed by higher order Lagrangians — namely by Lagrangians involving time 
derivatives of the degrees of freedom of order at least two — is an issue which 
by far has not been developed to the same extent[] as for systems whose 
Lagrangian only depends on the coordinates and their velocities^]. Never- 
theless, there do exist systems of physical interest described by such higher 
order Lagrangians, the most popular examples being perhaps higher order 
regularisations of quantum gauge field theories and so-called rigid strings [|], § 
or rigid particles^]. In fact, these examples involve the additional compli- 
cation that they possess local symmetries, leading therefore to constraints 
generating these gauge invariances on phase space. 

As is well known, there does exist a generalisation of the ordinary Hamil- 
tonian formulation in the case of higher order Lagrangians, which is due 
to OstrogradskyJF]. However, on the one hand, Ostrogradsky's approach is 
implicitly restricted to non constrained systems — which in particular do not 
possess local gauge invariances — , thus rendering this approach inapplicable 
to most, if not all physical systems of present fundamental interest. On the 
other hand, in Ostrogradsky's construction the structure of phase space and 
in particular of its local symplectic geometry is not immediately transparent, 
an obvious source of possible confusion when considering canonical or path 
integral quantisations of such systems. 

This note discusses how both problems can be resolved within the well es- 
tablished context of constrained systems^ described by Lagrangians depend- 
ing on coordinates and velocities only. Well known and powerful techniques 
become then immediately available, rendering the necessity of a separate dis- 
cussion of the quantisation of higher order systems — including constrained 
ones, and thus in particular their BRST quantisation — void of any justifica- 
tion. Any higher order system can always be cast in the form of an ordinary 
constrained system, namely one whose Lagrangian is a function only of first 
order time derivatives of the degrees of freedom, but not of time derivatives 
of higher order. 

That such a reduction of higher order Lagrangians is possible was indi- 

lr The classical analysis of such higher order regular or singular systems is available to 
some extent in the recent literature pjl. 
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cated already previously || |9|]. As should be clear, it suffices for this purpose 
to introduce auxiliary degrees of freedom associated to each of the succes- 
sive time derivatives of the original coordinates of the system. In effect, the 
canonical quantisation of rigid particles has already used|| the same idea, 
thus in a situation where Ostrogradsky's approach is not applicable as such. 

The present note is organised as follows. In the next section, Ostrograd- 
sky's construction is briefly considered. Sect. 3 describes how any higher order 
Lagrangian system can be cast into the form of a constrained system whose 
Lagrangian involves only first order time derivatives of the degrees of free- 
dom. The canonical Hamiltonian description of the auxiliary system is then 
addressed in Sect. 4 while its equivalence with Ostrogradsky's formulation is 
established in Sect. 5. Further comments are presented in the Conclusion. 

2 Ostrogradsky's Construction 

Let us consider a system with degrees of freedom x n {t) {n = 1,2, •••), t 
being the time evolution parameter of the system. Although the present 
analysis assumes that these coordinates are commuting variables, and that 
the index n takes a finite or an infinite number of discrete values, it should 
be clear that exactly the same considerations and the same conclusions as 
those developed hereafter are applicable to commuting and anticommuting 
degrees of freedom, as well as to an infinite non countable set of coordinates. 
The former case is that of bosonic and fermionic types of degrees of freedom, 
and the latter typically that of field theories. All the conclusions established 
in the present note are thus valid in complete generality, for any system 
described by some higher order Lagrangian. The restriction to a discrete 
set of commuting degrees of freedom is only one of ease of presentation. 
Moreover, in a first reading of the paper it might be useful to consider the 
case of only one degree of freedom||, namely ignore the index n altogether. 

These remarks having been made, let us assume that the dynamical time 
evolution of the system is determined from the variational principle being ap- 
plied to the action functional associated to some time independent Lagrange 
function 

Here, (m n > 1) (n = 1, 2, • • •) is the maximal order of all time derivatives of 
the coordinate x n (n = 1, 2, ■ ■ ■) appearing in the Lagrangian. In particular, 
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the discussion of this paper includes the familiar case when (m n = 1) for 
all degrees of freedom. Throughout the analysis, it might be interesting to 
consider the special case (m n — 1) (n — 1,2, •••) to see how well-known 
results are recovered from the present general situation. 

Note that the Lagrange function is assumed to depend on at least the first 
order time derivative of each degree of freedom x n . Otherwise, one would 
have to deal with some of the equations of motion being actually constraints, 
a situation not considered by Ostrogradsky. Moreover, the restriction to time 
independent Lagrange functions is again for reasons of convenience rather 
than of principle. Time dependent higher order Lagrangians can also be 
analysed along the same lines as developed hereafter. 

Considering the variational principle, it is clear that the Euler-Lagrange 
equations of motion of the system are given by 



n = l,2,.-.. (2) 



Following Ostrogradsky's lead[^] and in order to simplify the expression of 
these equations, let us introduce quantities p n ,a n ( a n = 0, l,---,m n — 1) 
defined recursively by 

Pn,i n — 1 = ~ 7~~T ~T,Vn,i n j = 1) 2, • • • , VTL n 1 , (3) 

OXn at 




with the initial value 



Pn,mn-\ — „ ( mn ) • (4) 



dXn 



The Euler-Lagrange equations of motion in (0) take then the simpler compact 
form 

which are very suggestive of Hamiltonian equations of motion. Note how 
these expressions generalise the familiar standard definitions in the case when 
(m n = 1) for all degrees of freedom x n (n = 1, 2, • • •). 

In order to reveal a possible Hamitonian description in the general case 
when the integers m n take arbitrary finite values, it is useful to consider 
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the differential of the Lagrange function Lq, in which the definitions of the 
variables p n ,a n are included. A little calculation then leads to the identity, 



m n — 1 



n o„=0 



n a „=0 



dL . 

"n Pi,0 



Note that the last sum in the r.h.s. of this expression is a combination of the 
Euler-Lagrange equations of motion of the system. 

The meaning of this result is as follows. Consider the quantity defined by 

m n — l 

h = E E 4 a "W« - Lo , (7) 

n a n =0 

with the variables p n ,a n determined by the recursion relations in (|3|) and . 
Since these variables involve the coordinates x n and their time derivatives 
up to a certain order which is different for each of the coordinates x n and 
each of the variables p n ,a n , so would a priori the quantity H given in (0). 
However, the identity (|]) establishes that this dependence is in fact rather 
specific, namely only through a dependence of the variables and p n ,a n 
(a n = 0, 1, • ■ • , m n — 1) themselves, 

H (x£ n) ,p n;an ) , a n = 0,l,---,TO„-l . (8) 

Note that this conclusion is valid irrespective of whether the Lagrangian 
Lq {x n , x n , • ■ • , x^" 1 "^ leads to constraints or not. 

The identity (|) also shows that the Euler-Lagrange equations of motion 
(H) are equivalent to the set of equations 

dH dH 

% n n = i Pn,a n = ~ 77~T > Oi n = 0, 1, • • • , Tfl n 1 . (9) 

9pn,a n dx { n n > 

In other words, the system of higher order Lagrangian Lq has been cast in 
Hamiltonian form, with the variables (a^™-*, Pn,m n ^j being canonically conju- 
gate pairs. Let us thus recapitulate. 
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Given the Lagrange function L (x n , x n , ■ ■ ■ , x^^j, first one introduces 
the conjugate momenta p n ,m n -i defined by 

dL 

Pn,m n -1 — ( mjj ) (%m %ni ' ' ' i X L "j 1 (^) 

Regarding the variables (x n , x n , ■ ■ ■ , x^ m " -1 ^ as independent, the relations 



fllOD may be inverted to give 

T (m») _ ~K-i) _ ™K-i) / T ~ ... T K-i) ^ ^ ni^ 

x n \-^ni^ni t-^n ! l J n,m n J ■ \- L - L J 

The remaining conjugate momenta jOn^-i (^n = 1> 2, ■ ■ • , m n — 1) are then 
determined by the recursion relations in (|3|). However, these relations are not 
used in order to express the variables x)^ n ~ 1 ' {i n = 1, 2, • ■ ■ , m n — 1) in terms 
of (x n , x n , ■ ■ ■ , a4 in_1 )), the conjugate momenta (p n ,i n i • ' ' ,Pn,m n -i) an d time 
derivatives of the latter. Indeed in Ostrogradsky's construction, the variables 
(x n , x n) ■ ■ ■ , a^" -1 )) have to be considered as being independent. We shall 
come back to this point shortly. 

Once the expressions for the quantities i^ 71 ™ -1 ^ determined as in (|TT|), the 
canonical Hamiltonian of the system is defined by (0), or equivalently, 



H(x£ n) 

„(m n -l) 



Pn,a„) 

m„-l 



,., 'Pn,m n -1 Lq (x n , X n , • • • , M + ^p n ^ n — i , (12) 



n i n =l 



in which only the substitutions for the variables x^ n ~ 1 ^ are performed, while 
the degrees of freedom x^ n ^ and p n ,a n («n = 0, 1, • • • , m n — 1) are considered 
as being independent. 

Finally, the Hamiltonian equations of motion are given by (|9|). Introduc- 
ing the fundamental Poisson brackets, 

{ X n ^ ' X m _ ^ ' { X n \Prn,am \ ~ ^nm^a n a m i {Pn,a n , Pm,a m } — , (13) 

with (a n = 0, 1, • • • , m n — 1) and (a m = 0, 1, • • • , m m — 1) (n,m = 1, 2, • • •), 
the equations ([]) take the Hamiltonian form 

4 Qn) = { x n n \ H } > Pn,a„ = {p n ,a n ,#} , « n = 0, m n - 1 . (14) 
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In other words, the variables (x^ n \p nja ^) (a n = 0, 1, • • • , m n — 1) are pairs 
of conjugate degrees of freedom, thus defining the phase space of the system 
and its local symplectic structure. 

Obviously, certain comments are in order. It is clear that Ostrogradsky's 
construction is applicable only to those higher order Lagrangians for which 
the inversions required in the determination of the quantities x^ mn_1 ^ are 
non degenerate. Namely, the Lagrangian L cannot lead to constraints of 
any kind. Constrained higher order Lagrangians are beyond the scope of 
Ostrogradsky's approach. 

Another issue with the present construction is the risk of confusion which 
arises when dealing with the variables x^ 1 ^ and their first order time deriva- 
tives i4 a ™) (a n — 0, 1, • • • , m n — 1), a situation which becomes even the more 
acute when considering canonical or path integral quantisations of such sys- 
tems. As emphasized above, only the first order time derivatives of the 
variables x^" 1 ™ -1 - 1 are to be solved for in terms of the conjugate momenta 
Vn,m n -i and the variables (a n = 0, 1, • ■ ■ , m n — 1), the latter considered 
to be independent of one another rather than being simply time derivatives 
of order a n of the coordinates x n . It is in this manner only that the canonical 
Hamiltonian defined in (12) can be made a function of the pairs of conjugate 
degrees of freedom (xi an s~Pn,«„)- To illustrate the possible confusion which 
might arise when this point is not fully appreciated, the reader is invited to 
consider a simple example in the case of a single degree of freedom x(t), such 
as, 

L Q (x,x,x) = -axx 2 — -bxx 2 , (15) 

with a and b being arbitrary constant parameters. If one attempts solving 
both for x and for x in terms of x and po and pi, there appear in the canon- 
ical Hamiltonian time derivative terms of the conjugate momentum p\\ It is 
thus important to develop Ostrogradsky's construction precisely in the man- 
ner emphasized above, keeping the variables x and x as independent, and 
inverting only for x in terms of x, x and p\. 

Nevertheless, when solving the Hamiltonian equations of motion ( |14D for 
the degrees of freedom x n (t), it becomes necessary, after having computed the 
Poisson brackets, to impose the condition that the variables x%™> are time 
derivatives of order i n (i n — 1, 2, • • • , m n — 1) of the coordinates x n (t). 

It is clear that both issues are solved at once by emphasizing explicitly 
the fact that in the Hamiltonian approach — hence also when considering 
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canonical and path integral quantisations of such systems — , all variables 
x (an) are kg considered as being independent. This is readily achieved by 
introducing independent auxiliary degrees of freedom, each corresponding to 
a time derivative of a given order of one of the original degrees of freedom. 
The same system can then be described in terms of an extended Lagrangian 
including a dependence on the auxiliary degrees of freedom, such that time 
derivatives of first order only are involved. In this manner, one is brought 
back||, § the realm of the usual type of dynamical systems for which 
most powerful techniques are available, with the additional advantage that 
constrained higher order Lagrangians do not need to be considered on a 
separate basis any longer. 



3 The Auxiliary Lagrangian 

Given a system of degrees of freedom x n (t) (n = 1,2, •••) with Lagrange 
function Lo(x n , x n , • • ■ , x^ 1 ^) (m n > 1) — be it regular or not — , let us intro- 
duce new independent variables q n ,a n (t) {ot n = 0, 1, • ■ ■ , m n — 1) such that the 
following recursion relations would hold, 

q n ,i n = q n ,i n -i , in = 1,2, • • • ,m n - 1 , (16) 
with the initial value 

<7n,0 = %n • (17) 

Clearly, the variables q n ^ n (i n = 1, 2, ■ • • , m n — 1) would then correspond to 
the time derivatives x^ of order i n of the coordinates x n , the latter being 
identical to the coordinates q n fl. 

In order to inforce the relations ( |I6| ) and (|T7j ) for the independent variables 
q n ,a n , additional Lagrange multipliers /i nin (t) {i n = 1,2, ■ • • ,m n — 1) are 
introduced. The variables (g n ,a n) /ViJ thus determine the set of independent 
degrees of freedom of the extended Lagrangian system, with (q n ,i n , Hn,i n ) 
{i n = 1,2,- ■ - ,m n — 1) being auxiliary degrees of freedom as compared to 
the original coordinates (x n (t) = q n fi(t)). The auxiliary Lagrange function of 
this extended description of the system is given by 

in) = 

mn— 1 

= A)(<?n,Cb qn,l, ■ ■ ■ , q n ,m„-l, 9n,m„-l) + X! X! (ln,i n ~ Qn,i n -l) ^n,i n ■ (18) 

n i n =l 
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Note that as advertised, the auxiliary Lagrangian L involves only first order 
time derivatives of the extended set of degrees of freedom. Obviously, due 
to the presence of the Lagrange multipliers [i n ,i n , the Lagrange function L 
in (0) defines a constrained system, to which the usual analysis |L0|, ||] of 



constrained dynamics is applicable. 

Before turning to that important issue however, let us first establish the 
equivalence of the auxiliary Lagrangian with the original formulation of the 
system determined by the Lagrangian L . Applied to L in ([18]), the varia- 
tional principle leads to the following equations of motion for the Lagrange 
multipliers n n)in , 

Qn,i n Qn,i n — 1 i ^"n lj 2, * • • , Tfl n 1 , (1^) 

while for the degrees of freedom q nyin (i n — 1, 2, • • • , m n — 1), one obtains, 



and 



9L d dL 

OC Ln,m n -l ™%,m„-l 

Finally, the equations of motion for q n are 

TT^- + An,l = . 22 

oq n ,o 

Note that the latter equations are in fact the actual equations of motion 
of the system. Indeed, all the other equations for q n ^ in and \i n ,i n are constraint 
equations which determine the auxiliary degrees of freedom q n ,i n in terms of 
successive time derivatives of the original coordinates (q n>0 = x n ), as well as 
the Lagrange multipliers /i n ,j n in terms of successive partial derivatives of the 
Lagrange function L . By substitution in ( p2]) of the successive definitions 
of the Lagrange multipliers ^ n ,i n i the equations of motion for g„ i0 reduce to 

sW^r^+t-wir/^o. , 23 , 

a n =0 \ dt J d( ln, an \dtj dq n . :mn -i 

Upon the substitution of the recursion relations (|19|) , one then indeed recovers 
the original Euler-Lagrange equations of motion in @). Hence, the complete 
equivalence between the auxiliary formulation of the system and the original 
one based on the higher order Lagrange function L (x n , x n , ■ ■ ■ , x^"^ is 
established. 
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4 The Hamiltonian Formulation 



Given the auxiliary Lagrangian formulation of higher order systems of the 
previous section, let us apply to it the ordinary analysis^] of constraints 
in order to develop its Hamiltonian description. The momenta canonically 
conjugate to the degrees of freedom q n , a „ (cn n = 0, 1, • • • , m n — 1) and ii n ,% n 
(i n = 1, , 2 • • • , m n — 1) are of course defined by, respectively, 

dL dL 

Pn,a n — j ^n,i n — T7T ■ (24) 

u Hn,a n < JfJ'n,i rl 

However, the phase space degrees of freedom (q n ,a n iPn,a„] /V«™> ^n,in) are 
not all independent. In fact, the system possesses the following primary 
constraints, 

$n,i„ = , Tt n ,i n = , i n = 1, 2, • • • , m n - 1 , (25) 

where 

$n,in = Pn,i„-1 + fJ>n,i n > *n = 1, 2, • • • , m n - 1 . (26) 

Both sets of primary constraints follow from the particular way in which the 
auxiliary degrees of freedom are introduced in the definition of the extended 
Lagrange function L in flT8|). The primary constraints obey the algebra of 
Poisson brackets 

{^n,in> ^m,j m } ' {^n,i n > ^m,i m } ^nm^i n i m i \^n,in ] ^m,im} , (27) 

with (i n = 1, 2, • • • , m n - 1) and (i m = 1, 2, • • • , m m - 1) (n,m = 1, 2, • • •), 
showing therefore already at this stage that the primary constraints are cer- 
tainly also second class constraints. 

Among all conjugate momenta, p n ,m n -i certainly play a distinguished role 
since on the one hand, they are the only ones not involved in any of the pri- 
mary constraints above, and on the other hand, their conjugate coordinates 
?n,m„-i are the only variables whose first order time derivatives do appear in 
the original Lagrange function L . Indeed, we have 

^Lq / ' \ / \ 

Pn,m„-1 = ~7p l?n,0> Qn,lj * ' ' j Qn,m n -1, Qn,m„-l) ■ \^°) 
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As in Ostrogradsky's approach, let us then assume that for fixed values of 
<?n,«„ (a = 0, 1, • • • , m n — 1), these relations are invertible, leading therefore 
to the velocities, 

Qn,m n — 1 = Qn,m n — 1 {Qn,a n > Pn, m n — l) • (29) 

In other words, given fixed values for q n ,i n -i (i n = 1,2,- ■ ■ ,m n — 1), the 
dynamical system of degrees of freedom q n ,m n -i with Lagrange function 
Lq (q n ,o, q n ,x, ■ ■ ■ , Qn,m n -i, q n ,m n -i) is assumed to be a regular system, namely 
not leading to any constraints for the conjugate momenta p n ,m„-i- Conse- 
quently, the constraints in determine the full set of primary contraints in 
the extended formalism of the higher order Lagrangian L ( 

The distinguished role of the conjugate variables (q n ,m n -i, Pn,m n -i) justi- 
fies the definition of the restricted Legendre transform of L (q n ,a n , ?n,m„-i)> 
leading to the restricted canonical Hamiltonian, 

Ho(q n ,a n iPn,m n -l) = 2j Qn,m n -lPn,m n -l ~ L (q n ,a n , ?n,m„-l) • (30) 

n 

In the same way as was established for the Hamiltonian H in (^), note that 
the restricted Hamiltonian H is a function of the variables (q n ,a n ,Pn,m n -i) 
only, irrespective of whether the relations are invertible or not, namely 
irrespective of whether L (q n>0) q n>1 , ■ ■ ■ , q n , mn -u ?n,m„-i) defines a regular 
system in the coordinates q n ,m n -\ or not0. In the present discussion, the 
assumption of regularity is necessary only in order that no further primary 
constraints beyond those in ([25|) appear in the analysis. 

In terms of the definitions and the primary constraints above, the cano- 
nical Hamiltonian of the extended system, 

m„ — 1 m n — 1 

H = J2 £ Qn,a n Pn,a n + E E ^n,i n ^n,i n ~ L , (31) 
n a n n i n =l 

is readily found to be given by 

m n — 1 

-^0 {qn,a n ,Pn,a n 'i fJ>n,i n ) = #0 {<ln,a n , Pn,m n -l) ~ X! ^n»*n9n,i„- (32) 

n j n =l 

However, as is well known ||, due to the presence of constraints, the Hamil- 
tonian generating the genuine time evolution of the system under which the 



10 



constraints are preserved, is in general given by the canonical Hamiltonian 
Ho and a linear combination of the constraints. Thus in the present case, 
the would-be Hamiltonian is of the form, 



m n — 1 



n i„=l 



(33) 



with A^j n and A^ n being Lagrange multipliers for the constraints. Consistent 
time evolution of the primary constraints $ n ,i n and n n ,i n then imposes the 
relations 



A 



(i) 



1, 2, • • • , Tfl n 1 , 



(34) 



as well as 



(2) <9# ,(2) 
<9<?n,0 



- , in = 2, 3, ■ ■ • , m n - 1 . (35) 



Consequently, the extended formulation of the system does not possess se- 
condary constraints, while its extended Hamiltonian reduces to 

mn — 1 

F* = Hq + £ Qn,i„Pn,i n -l + 

n i n =l 



+ E^ 



dH 

dq n ,c 



m.n — 1 



E E 



7T 



dH r , 



dq n 



(36) 



n j„=2 

However, as already pointed out previously, all primary constraints $ n ,i n 
and Tc n ,i n are second class, and may thus be solved for explicitly provided the 
canonical Poisson brackets are traded for appropriate Dirac brackets ||T0|, §]. 
Choosing to solve the constraints in terms of 



r"n,if, 



, • • • , m n - 1 



(37) 



the reduced phase space degrees of freedom are then simply (q n ,a n ,Pn,a n ) 
(a n = 0, 1, ■ • ■ , m n — 1). On the other hand, given the algebra (|27J) of con- 
straints, the Dirac brackets of the reduced Hamiltonian description are easily 
seen to remain canonical (a n = 0, 1, ■ • • , m n — 1; a m = 0, 1, ■ • • , m m — 1), 

{?»i,a„i ?m,Q m }fl , {Qn,a n i Prn,ctm} D ^nm^a n a m j {Pn,a n > Prn,a m } D ■ 

(38) 
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The constrained description of higher order Lagrangian systems has thus 
lead to the following Hamiltonian formulation. The local phase space coordi- 
nates are the canonically conjugate pairs of degrees of freedom (q n ,a n ,Pn,a„) 
(a n = 0, 1, • • • , m n — 1), with the local symplectic structure determined by 
the Dirac brackets in (|3^). Time evolution in phase space is specified through 
these brackets by the extended Hamiltonian, 

m„— 1 

He (q n ,a n ,Pn,a n ) = H Q {q n ,a n , Pn,m n -l) + Qn,inPn,i n -l ■ (39) 

n i n =l 

In particular, the fundamental equations of motion are (i n = 1, 2, • • • , m n — 1), 

Qn,in—X Qn,i n > <?n,m n -l ~n j l^UJ 

OPn,m n — l 

as well as 

dH . dH 

Pn,o = —5 ' 3V<» = ~a Pn,i n -1 • (41) 

However, the restricted canonical Hamiltonian Hq is such that 

,m n - l) = ~a (in q„ Pn ,m n — 1 )) , (42) 



so that the equations fl4"I|) for p n ,a„ are equivalent to 

" = ■ (43) 

with 

Pn,i n -i = ^— — - -r,Pn,i n , z„ = l,2,---,m n - 1 . (44) 

%M„ at 

Expressed in this manner, it is clear how these Hamiltonian equations 
of motion are indeed equivalent to the Euler-Lagrange equations (|i~9|) to 
( |22|) of Sect. 3, when the constraints (/i n ,i„ = — Pn,i„-i) are accounted for. 
Indeed, ([43"D determine the actual equations of motion of the system and 
are equivalent to (f22"D. The conjugate momenta p Ut0 are defined recursively 
through the equations fl44|), which are equivalent to ( p0|) and ([211), given the 
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momenta p n mn _i. Finally, the latter quantities, which are absent of course 
in the Lagrangian equations of motion except through their definition as 

Vn.m n -\ = 777 (Qn,a n i Qn,m n -l) > (^5) 

OQn,m n — l 

are determined implicitly by the Hamiltonian equations of motion, 

Qn,m n — 1 = 77 (?n,a n i Vn,m n — l) • (46) 

Since the Lagrangian L Q (q n ,a n , Qn,m n -i) is assumed to be regular in the co- 
ordinates g n)mn _i, this latter relation is indeed invertible, leading back to 
the relation (|45|). It is in this way that the present Hamiltonian equations of 
motion are equivalent to the Euler-Lagrange equations under the Lagrangian 
reduction, namely the reduction of conjugate momenta p n ,a„ in terms of the 
coordinates g n Q , n and their velocities g n ,a„- Therefore, since the auxiliary 
Lagrangian formulation was shown to reproduce the Euler-Lagrange equa- 
tions of the higher order Lagrangian, the present Hamiltonian construction 
is established to be equivalent to the original description of the system as 
well. 



5 Ostrogradsky's Approach Revisited 

The equivalence of the results obtained through the analysis of constraints 
applied to the auxiliary formulation of regular higher order Lagrangian sys- 
tems with Ostrogradsky's construction is now obvious. 

In the latter approach at the Hamiltonian level, the successive time deriva- 
tives («n = 0, 1, • ■ • , m n — 1) of the degrees of freedom x n have to be con- 
sidered as being independent. In the constrained formulation, these variables 
correspond to the independent auxiliary coordinates q n ,a n , with in particular 
(<7n,o — x n)- In addition, the fundamental brackets ( |13"D in Ostrogradsky's 
formulation are identical to the canonical Dirac bracket s(|3~8|) of the reduced 
phase space degrees of freedom (q n ,a n ,Pn,a n )- 

Finally, it is clear that the extended Hamiltonian He in ( |39D is identical 
to Ostrogradsky's canonical Hamiltonian H in (|T2|). In particular, note how 
in the definition of the latter quantity, the restricted canonical Hamiltonian 
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H defined in ( |30| ) appears naturally, indeed emphasizing once again the 
distinguished role played by the time derivatives x^""^ of maximal order of 
the degrees of freedom x n (t). 

Therefore, the analysis of the previous section, based on the auxiliary for- 
mulation of regular higher order Lagrangian systems and Dirac's analysis of 
constraints, has recovered precisely Ostrogradsky's Hamiltonian description 
of such systems. 

6 Conclusion 

This note has established the equivalence of the Hamiltonian formulation of 
regular higher order Lagrangian systems due to Ostrogradsky0, with a con- 
strained auxiliary description^], £5] of such systems in which time derivatives 
of degrees of freedom of at most first order only are involved. The latter 
approach offers the following advantages, however. 

In Ostrogradsky's construction, time derivatives of the coordinates of dif- 
ferent order have to be considered as being independent. Such a situation is 
a possible source of confusion, especially at the quantum level when translat- 
ing Poisson brackets into ( ant i) commutation relations for the fundamental 
quantum operators. Indeed, it is not always clear when to consider a time 
derivative of given order as an independent variable or as the first order time 
derivative of some other variable in Ostrogradsky's phase space. In the aux- 
iliary approach, this issue is avoided altogether ab initio, since independent 
auxiliary degrees of freedom are introduced explicitly, each being associated 
with a time derivative of given order of each of the original degrees of freedom. 
In this manner, the local structure of phase space and its local symplectic 
geometry, is made perhaps much more transparent than in Ostrogradsky's 
approach. 

More importantly however, the auxiliary formulation presents the addi- 
tional advantage that the auxiliary Lagrangian depends on time derivatives 
of first order only. Therefore, any higher order Lagrangian system — be it 
regular or not — can always be brought into the realm of those Lagrangian 
systems for which a wealth of methods — classical and quantum — have been 
developed over the years. Due to the presence of auxiliary degrees of free- 
dom, the auxiliary formulation always leads to constraints, requiring the 
techniques of constrained dynamics 0]. 
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Finally, in contradistinction to Ostrogradsky's construction which ap- 
plies to regular higher order systems only, the auxiliary formulation, being 
already a constrained one, does not require to distinguish between regular 
and singular higher order Lagrangian systems. Hence, the quantisation of 
such systems, including the BRST quantisation of singular ones, does not 
necessitate a separate and generalised formalism not yet developed. All the 
readily available methods of ordinary constrained quantisation — and nothing 
more — suffice for the Hamiltonian formulation and the quantisation of any 
higher order Lagrangian system. As this note has established, Ostrogradsky's 
construction is thereby recovered exactly in the case of regular systems. The 
case of singular systems however, is beyond the scope of the latter approach, 
and the auxiliary formulation then becomes unavoidable. In effect, precisely 
this method has been applied already to rigid particles for example, with 
important conclusions as to their quantum consistency [[]]. 
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